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Abstract 

Despite decades of progress, an understanding of unconventional superconductivity still remains 
elusive. An important open question is about the material dependence of the superconducting prop¬ 
erties. Using the quasiparticle self-consistent GW method, we re-examine the electronic structure 
of copper oxide high-Tc materials. We show that QSGW captures several important features, dis¬ 
tinctive from the conventional LDA results. The energy level splitting between d„, 2 _y 2 and d^^ 2_^2 
is significantly enlarged and the van Hove singularity point is lowered. The calculated results com¬ 
pare better than LDA with recent experimental results from resonant inelastic xray scattering and 
angle resolved photoemission experiments. This agreement with the experiments supports the pre¬ 
viously suggested two-band theory for the material dependence of the superconducting transition 
temperature, Tc. 
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After the seminal work of finding high temperature superconductivity in a ceramic copper 


Q. 


oxide material llH, efforts to understand the cuprate have long been a central part of mod¬ 
ern condensed matter physics. Although many intriguing aspects of its electronic behaviors 
have been unveiled and great progress has been made, an understanding of its supercon¬ 
ducting mechanism, the novel interplay between competing phases, and its relationship to 
other correlated phenomena is still far from clear {2 12|. One simple, well-dehned but still 
open question is what determines the superconducting transition temperature (Tc), or its 
material dependency. For example, the Tc of the single layer cuprates can be different by 
a factor of two; ~40K for La 2 Cu 04 and ~90K for HgBa 2 Cu 04 . On the one hand, it may 
be too early to ask this question, while the superconducting mechanism itself still remains 
elusive. On the other, however, figuring out the detailed features behind the material de¬ 
pendency can provide the crucial hint for further understanding of superconductivity and 
other related properties. In fa^ many directly or indirectly related theoretical studies have 
been performed on this issue 13l-l25|. 

In this regard, a notable suggestion has recently been made 26l-l29| . According to this 
theory, the Tc of the single-layer cuprate can be described with a two-orbital model that 
considers both the and Wannier orbitals, and the energy level offset between 

the two orbitals, AE, plays a key role in determining Tc. Whereas the larger value of this 
energy separation produces the higher Tc {e.g., the case of HgBa 2 Cu 04 ) due to the better one 
band feature achieved, the smaller value results in the lower Tc {e.g., the case of La 2 Cu 04 ) in 
spite of the better nested Fermi surface. The calculated Eliashberg parameter (A) based on 
the many-body calculation using the fluctuating exchange (FLEX) approximation 30|, 1M| 
clearly exhibits a linear dependence on AE while the other parameters are shown to be less 


important. Further, this theory can be extended to t 
correlation between the Fermi surface shape and Tc 


re bilayer case 


29| . which explains the 


24|. 


A possible experimental test to verify this two-band theory is to examine the correlation 
between Tc and AE, the Fermi surface shape, or the partial density of states of the 
orbital, which can be measured by recent techniques such as resonant inelastic xray scattering 
(RIXS), and (angle resolved) photoemission spectroscopy (ARPES). However, while the 
theoretical AE or the Fermi surface shape was obtained from the LDA and used as the 


‘inputs” for the many-body FLEX calculation in Refs. 


26 


29 I. the experimentally determined 


AE and the Fermi surface shape should be regarded as the “outputs” or “results” after the 
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consideration of the many-body correlation effects beyond LDA/GGA. In fact, while a RIXS 
study reports that Tc is higher for larger AE 3^, the actual experimental value of AE is 
larger than the theoretical evaluation, presumably due to this “input vs. output” problem. 
One possible way to resolve this problem, at least partially, is to re-evaluate AE as an output 
of the FLEX calculation. However, this approach would suffer from various ambiguities 
regarding the Hubbard interaction strength and the dehnition of the renormalized AE. 
It is problematic since a quantitative comparison is required in between the th eory and 
experiment, while only the qualitative comparison was made regarding Tc in Refs. 


26 


29 


In the present paper, we use a hrst-principles approach, exploiting the quasiparticle self- 
consistent GW (QSGW) method. It enables us to take into account the correlation effects 
beyond LDA/GGA. In this way, we can obtain a well-dehned renormalized AE without 
introducing adjustable parameters. 

In Ref. [iTl, a quantum chemical approach was adopted to evaluate the energy level offset 
between Cu-dx 2 -y 2 and orbitals, where the correlation effects were taken into account 

within a cluster-based conhguration-interaction-type calculation. A good agreement with 
the RIXS experiment was found by assuming the energy difference of ferromagnetically and 
antiferromagnetically ordered states to be 2J, where J is the antiferromagnetic coupling 
constant. Our approach is fairly different and is along the line of the hrst-principles band 


calculation as in Refs. mm 


. £- 


In the sense mentioned above, the calculated AE can be 


compared to the experiments, while it should not be regarded as an input parameter for 
the many-body calculation, because doing so would result in a partial double counting 
of the correlation effects. Still, the present approach can also provide a hrst-step hint 
toward obtaining a better “non-interacting” Hamiltonian that can be used as an input 
for the many-body calculation of superconductivity. In fact, it is known that the non¬ 
interacting Hamiltonian obtained from LDA has a problem when used as an input for the 
FLEX calculation, and the LDA/GGA estimation of AE for La 2 Gu 04 is too small to account 
for the maximum Tc of 40K in the La 2 Gu 04 33(]. In this context, it is worth pointing out that 


the GW method has been successfully applied to the many of strongly corre 


in combination with, for example, dynamical mean held theory (DMFT) 


34 


ated materials 


37|. 
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Results and Discussion 


To our knowledge, there is no previous QSGW study for the cuprate band structure 
although it has been discussed conceptually |38| . Here we first examine the electronic struc¬ 
ture and the two-band theory for the material-dependent of a single layer cuprate. While 
the QSGW calculation produces notable differences in the band structure and Fermi surface 
from LDA, the two-band explanation for still remains valid. QSGW results of model 
parameters are presented and compared to the RIXS data as well as the LDA calculations. 
It clearly shows that the parameters produced by QSGW are in better agreement with the 
experiment. Finally, we investigate the epitaxially strained La 2 Gu 04 whose noticeable 
increase has been previously reported. Two-band theory also works well for this situation. 


Electronic structure and the Tc of single layer compounds 


Fig. [U^a) and (b) show the band dispersion and projected density of states (PDOS) of 
La 2 Gu 04 calculated by LDA and QSGW, respectively. The LDA result is in good agreement 


with the previous calculation (see, for example. Ref. 


26 


and Ref. 


271 ). In the QSGW, several 


important differences are noted. First, the band width of both Cg orbitals are significantly 
reduced, by about 1.30 and 0.65 eV for dx 2 -y 2 and respectively, indicating that the 

band width overestimation (or effective mass underestimation) problem of LDA is somehow 
overcome by the QSGW procedure. Another key difference is that the separation between 
the dx 2 -y 2 and bands becomes larger in QSGW, as seen in Fig. [11(b). The d^^-y^ 

energy level, Ea, 2 _y 2 , is shifted from —0.14 (LDA) to —0.03 (QSGW) while 'Ej^z 2_^2 is from 
—0.98 (LDA) to —1.68 (QSGW) (indicated by the arrows in Fig. [H see also TableH]). It is a 
factor of two difference in AEe ; 0.84 eV in LDA and 1.66 in QSGW. The correct estimation 


of this quantity is important especially in the two-band theory for T,. 


26 


291. The large 


value of AEp is indeed consistent with the RIXS data as will be discussed further below. 

The same features are also found in HgBa 2 Gu 04 , as presented in Fig. [T](c) and (d). The 
dx 2 _y 2 band width is reduced by ~0.75 eV in QSGW compared to LDA and its center 
position of PDOS moves slightly upward by 0.15 eV. While the d^z^_r 2 dispersion in this 
material is already quite small due to the thicker blocking layer, its band width in QSGW 
is further reduced. AEg is 2.25 in QSGW, again noticeably larger than the LDA value of 
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1.67 eV (see Table [I]). 

In QSGW, the of both La 2 Cu 04 and HgBa 2 Cu 04 is enhanced compared to 

LDA/GGA. How does this affect the theoretical estimation of T^? First, it should be noted 
that these parameters cannot be directly adopted as the inputs for the FLEX calculation. 
This is because, in principle, the QSGW self energy should be partially subtracted before we 
put it into any of many-body calculations. While there is no well-dehned prescription yet 


for this kind of ‘double-counting’ problem 


, the “best” AEg that should be adopted in 


the FLEX evaluation of Tg may be lying somewhere in between the QSGW and LDA/GGA 
values. This can provide better quantitative agreement with the experiment, especially in 
La 2 Gu 04 , for which the LDA/GGA value of AEg^ is found to be too small to account for 
Tg = 40K. 

Oxygen states are also affected. Gompared to LDA results, the 0-2p levels obtained by 
QSGW are signihcantly lowered in energy, as indicated in Fig. |2l As summarized in Tabled 
the center position of in-plane oxygen PDOS is located at —4.44 (—3.55) eV in LDA and 
at —5.06 (—4.05) eV in QSGW for La 2 Gu 04 (HgBa 2 Gu 04 ). The same feature is found for 
the apical 0-pz PDOS. As a result, the energy difference, AEp=EapicamEinpiane, is changed 
from 1.68 (1.00) in LDA to 0.55 (—0.86) in QSGW for the case of La 2 Gu 04 (HgBa 2 Gu 04 ), 
see Table [T1 The correct estimation of AEp is also important for understanding Tg since it i 




IS 


an underlying quantity to determine AEg^ AE AE^ -|- AEp [^) in combination with 
other parameters. 

Some other changes produced by QSGW are also noted. The (iszS-gZ components in the 
bands below —1.5 eV in Fig. [11(a) are reduced in QSGW, and the free-electron-like bands 
at P and Z points above the Fermi energy are shifted upward. As the position of the t 2 g 
complex is lowered (red color), the band has almost no mixture with other bands 

below the Fermi energy. Higher-lying La-4/ bands (not shown) move further upward as has 


been previously noted in the nickelate systems 


kll- 


Fermi surface 

The shape of the Fermi surface is important for understanding cuprate superconductivity. 
For example, its nesting is crucial for the spin fluctuation pairing. Also, a notable correlation 
between the experimentally observed Tg at the optimal doping (T™’^’^) and the Fermi surface 
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warping has been identified by Pavarini et al. 2^. Here we discuss the Fermi surface 
calculated by QSGW in comparison to the LDA result and experiment. 

The calculated Fermi surfaces are presented in Fig. [21 LDA ((a, c)) and QSGW ((b, d)). 
The hole doping is simulated by the rigid band shift method so that the electron occupation 
in Cg orbitals is reduced by O.lSe per unit cell. Notable features are found in the QSGW 
Fermi surface for La 2 Gu 04 . Gontrary to the LDA result of Fig. [3]^a), Fig. [3]^b) has the pocket 
centered at (tt, tt) point as in HgBa 2 Gu 04 Fermi surface (see Fig. [3](c) and (d)). This feature 
is in good agreement with ARPES data ^ which also reports the pocket centered at (tt, 
tt) point. Further, the (i 3 ^ 2 _^ 2 -orbital character (dark purple) is significantly reduced and 
the dx 2 _y 2 character (bright yellow) is dominant in the QSGW result, which is distinctive 
from the LDA in which the significant amount of d^z 2_^2 components are observed near (vr, 
0) and (0, vr). 

In the case of HgBa 2 Gu 04 , the difference between LDA and QSGW is less pronounced, 
see Fig. n c) and (d). While the QSGW Fermi surface is slightly more rounded, the overall 
shape is not much different. Since the dx' 2 -y 2 orbital character is dominant and 
band is well separated from Fermi level already in LDA due to the thicker blocking layers 
enhancing two-dimensional feature, the LDA result is quite similar to QSGW. 


Comparison with RIXS 


We now turn to the comparison with the RIXS data. Recently, Sala et al. |32j successfully 
extracted the important model parameters for several different cuprate materials based on 
RIXS spectra. In this subsection, we examine the material dependent parameters by QSGW 


and compare them to the experimental values. With Ref. 


32 


two more compounds, namely, Sr2Gu02Gl2 and GaGu02 [43 1. 


as our main reference, we include 


Our results are summarized in Table lH Fig. IH and Fig. [5l The values of lODq, defined 
as the difference between two energy levels of dx‘ 2 -y 2 and d^y (see Fig. 111(a)), are larger in 
the QSGW calculation by ~62-125% than the LDA values. While LDA underestimates 
lODg compared to the experiment, QSGW slightly overestimates, which is related to the 
tendency that Cn-t 2 g bands are pushed down relative to e^, as was also observed in the 


previous QSGW calculations for other transition-metal oxides 


41 


44 


45| . It is important 


to note that overall the QSGW result is in better agreement with experiment, as clearly 
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seen in Fig. lU^b). 

As noted in the above, according to the two-band theory by Sakakibara et al. 


26 


m 


the important parameter that governs is (or ADs+bDt in Ref. l32l) . Fig. 111(b) and 
Fig. m clearly show that the calculated values of AEg^ by QSGW are in excellent agreement 
with those from RIXS spectra; the difference is 2-8 %. The LDA values are noticeably 
smaller than the experiments although the difference gets reduced in the higher Tc materials, 
CaCu02 and HgBa 2 Cu 04 (see Fig. [5]). This can be taken as a strong support for the two- 
band theory in the sense that the LDA value of AEe^ as an input for FLEX provides 
qualitative information of material dependence, while the AEe^ by QSGW already contains 
the correlation effect beyond LDA, being consistent with RIXS. 


Another parameter deduced from RIXS in Ref. 
tween dxy and dyz,zx- In this 
and experiment (see Fig. |4](b)). 


32l is 3Ds-5Dt, the energy level difference 


between dxy and dyz,zx- In this case, the LDA results are not much different from QSGW 


The effect of epitaxial strain 

An interesting aspect found in the Tc trend of the cuprates is its signihcant enhancement 
in the thin him form. Locquet et al. reported 46j that Tc can be controlled by epitaxial 


strain by about factor of two 


47( 1 . The underdoped La 2 Gu 04 with its bulk Tc of 25K exhibits 


a higher and lower Tc of ~49 K and 10 K when it is grown on SrLaA 104 (SLAG) and SrTiOs 
(STO) substrates, respectively |46|. It is therefore important to check whether the two-band 
theory is also consistent with this observation. 

In order to simulate the tensile and compressive strain produced by STO and SLAO, 
we hrst optimized the c lattice parameter with two different in-plane lattice constants, 
aSTO =3 9 Q 5 and a®^^°=3.755 A, for La 2 Gu 04 , which originally has ao=3.782 A and 
Co=13.25 A. As expected, the optimized out-of-plane parameters get smaller and larger 
under the tensile and compressive strain, respectively; CgTQ=12.96 and CgLAo = 13.36A. As 
a result, the ratio between the out-of-plane and in-plane Gu-0 distance, r = dapicai/Anpiane, 
is found to be 1.32, 1.28, and 1.24, for Oq and a®™, respectively. 

The calculated values of AEgg are plotted in Fig. [6l Both LDA and QSGW predict that 
AEgg gets enhanced and reduced under compressive and tensile strain, respectively, which 
is consistent with the experimental observation The reduction of AEg at a=asTO is 
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about 0.16 eV in both LDA and QSGW, and the enhancement at a=asLAo is 0.29 (LDA) 
and 0.47 eV (QSGW). 


Summary and Conclusion 


Using the QSGW method, we re-examined the electronic structure of copper oxide high 
temperature superconducting materials. Several important features were found to have been 
captured by the GW procedure, such as effective mass enhancement. The shape and orbital 
character of the Fermi surface were also notably changed, especially for the case of La 2 Gu 04 , 
and they are in good agreement with the ARPES data Important model parameters 
including the key quantity for the two-band theory of Tc, AEe^, were examined, and the 
QSGW results were in excellent agreement with RIXS data. 

The present study shows that the first-principles band calculation can quantitatively re¬ 
produce the experimental observation by taking into account the correlation effects beyond 
LDA. We emphasize that it is not inconsistent with the previous study by Sakakibara et 
al. which takes the LDA result as an input for the many-body calculation of supercon¬ 
ductivity. While the QSGW result cannot be used as a direct input for the FLEX-type 
calculation because of the partial double-counting of the many-body correlation, the “best” 
non-interacting Hamiltonian, that can serve as an input, may lie somewhere in between 
the LDA and QSGW. Obtaining a well-dehned non-interacting Hamiltonian is, therefore, 
an important future direction for the hrst-principles-based description of high-temperature 
superconductivity, and it may quantitatively resolve the problem of low Tc in La 2 Gu 04 
produced by the LDA input [SSj. 


Methods 


Quasiparticle self-consistent GW 


The QSGW l4a, l48| calculates Hq (non-interacting Hamiltonian describing quasipar¬ 


ticles or band structures) and W (dynamically-screened Goulomb interactions between the 
quasiparticles within the random phase approximation) in a self-consistent manner. While 
the ‘one-shot’ GW is a perturbative calculation starting from a given Hq (usually from 









LDA/GGA), QSGW is a self-consistent perturbation method that can determine the one- 
body Hamiltonian within itself. The GW approximation gives the one-particle effective 
Hamiltonian whose energy dependence comes from the self-energy term S(ci;) (here we omit 
index of space and spin for simplicity), and in QSGW, the static one-particle potential 
is generated as 


? E I*) + R<!|E(£j))„} (^,1, 


( 1 ) 


V 


where e* and refer to the eigenvalues anc 


is the Hermitian part of the self-energy 


44 


eigenfunctions of ffo, respectively, and Re[S(£:)] 


45 


48|. With this one can define a new 


static one-body Hamiltonian Hq, and continue to apply GW approximation until converged. 
In principle, the hnal result of QSGW does not depend on the initial conditions. Previous 
QSGW studies, ranging from semiconductor s [ 44 

oxides 4^, l45|, and 4/-electron systems [5^, have demonstrated its capability in the 


45| to the various Sd transition metal 


description of weakly and strongly correlated electron materials. 


Computation details 


We used our new implementation of QSGW [Ml by adoptin g th e ‘augmented plane wave 
(APW) -|- muffin-tin orbital (MTO)’, designated by ‘PMT’ 52, 53|, for the one-body solver. 
The accuracy of this full potential PMT method is proven to be satisfactory in the supercell 
calculations of homo-nuclear dimers from H 2 through Kr 2 with the significantly low APW 
energy cutoff of ~ 4 Ry, by including localized MTOs 53| . A key feature of this scheme for 


QSGW is that the expansion of can be made with MTOs, not APWs, which enables us 
to make the real space representation of at any k point. 


We performed the calculations with the experimental crystal structures 54j-l57l|. and 
used 10x10x10, 12x12x12, 12x12x8, and 14x14x14 k points for LDA calculations of 
Sr2Gu02Gl2, La 2 Cu 04 , HgBa 2 Cu 04 , and CaCu02, respectively. As for QSGW calcula¬ 
tions, in order to reduce the computation cost, the number of k points were reduced to be 
5x5x5, 7x7x7, 8x8x4, and 8x8x8 for the hrst Brillouin zone of Sr2Cu02Cl2, La 2 Cu 04 , 
HgBa 2 Gu 04 , and GaGu02, respectively. The MTO radii used in our calculations were as 
follows: (i) 1.58, 1.04, 0.89, and 1.38 A for Sr, Gu, O, and G1 in Sr2Cu02Cl2, (ii) 1.43, 0.97, 
and 0.86 A for La, Cu, and O in La 2 Cu 04 , (iii) 1.10, 1.59, 1.05, and 0.83 A for Hg, Ba, Cu, 
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and O in HgBa 2 Cn 04 , and (iv) 1.54, 1.01, and 0.86 for Ca, Cn, and O in CaCn02. 


Many of the key parameters in this study are dehned in terms of the energy levels of each 
orbital, such as 'Ejx2_y2 and 'Ej^^2_^2. To quantify them we simply take the center of mass 
position of PDOS: 


St:: Egc(E)dE 

Sl~g^(E)dE ' 

^ -^min 


( 2 ) 


where ga{E) is PDOS for a given orbital a. An ambiguity is inevitably introduced in 
determining Emin max, and we set the range to cover the whole antibonding band complex 
for Cn-Cg states. Emax) for La 2 Cu 04 is (—1.95 eV, 2.05 eV) in PDA and (—2.20, 1.55) 

in QSGW. For HgBa 2 Cu 04 , the band dispersion changes and the values of and Emax 
are redehned accordingly: (E((^;^, Emax) = (—2.40, 2.50) in PDA and (—2.55, 1.65) in QSGW. 


Importantly, none of the reasonably dehned energy ranges change our conclusion 


, and 


the values are well compared with those reported in the previous study using a maximally 


localized Wannier function 


26 


27|. 
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Tables 


TABLE I: The calculated parameters by LDA and QSGW. The px,y is from the in-plane oxygen, and 
Pz from the out-of-plane. The definitions of parameters can be found in Fig. HI The experimental 
values are taken from Ref. 32. 



Epermi (eV) Erf^j “ EPermi (eV) '^Op^ y “ 

Epermi (eV) “ 

Epermi (cV) 


LDA 

QSGW 

LDA 

QSGW 

LDA 

QSGW 

LDA 

QSGW 

Sr2Gu02Cl2 

-0.29 

0.01 

-1.23 

-2.05 

-4.67 

-4.92 

-3.61 

-5.44 

La 2 Cu 04 

-0.14 

-0.03 

-0.98 

-1.68 

-4.44 

-5.06 

-2.76 

-4.06 

HgBa 2 Gu 04 

-0.20 

-0.05 

-1.87 

-2.30 

-3.55 

-4.05 

-3.00 

-4.91 

GaGu02 

-0.30 

-0.17 

-2.34 

-2.63 

-3,67 

-4.05 

None 

None 


AEe„ (4Ds + 5Dt) (eV) 


3Ds — 5Dt (eV) 


lODq (eV) 


LDA QSGW 


Exp 


(Ref. 


3 ^ 


LDA QSGW 


Exp 


(Ref. 




LDA QSGW 


Exp 


(Ref. 


M) 


Sr2Gu02Cl2 

0.94 

2.05 

1.97 

0.22 

0.21 

0.33 

0.81 

1.82 

1.50 

La 2 Cu 04 

0.84 

1.66 

1.70 

0.17 

0.11 

0.32 

1.47 

2.13 

1.80 

HgBa 2 Gu 04 

1.67 

2.25 

None 

0.42 

0.42 

None 

0.90 

1.57 

None 

GaGu02 

2.04 

2.46 

2.65 

0.05 

0.28 

0.31 

0.92 

1.69 

1.64 
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Figure Legends 

FIG. [U The band dispersion and PDOS of La 2 Cu 04 (a, b) and HgBa 2 Cu 04 (c, d) 
calculated by (a, c) LDA and (b, d) QSGW. The green and blue colors refer to the dx 2 _y 2 
and characters while the size of the colored dots represents their weight. The same 

color scheme was used for PDOS. Other bands than the two eg states are represented by 
red color. The center of mass position of PDOS is marked by an arrow. Fermi energy is set 
to be 0. 


FIG. 121 The calculated oxygen PDOS of La 2 Gu 04 by LDA and QSGW. (a) In-plane 
oxygen p^ y and (b) out-of-plane oxygen orbitals are plotted. The center of mass position 
of each PDOS is marked by an arrow. Fermi energy is set to be 0. 


FIG. |3l The orbital-resolved Fermi surfaces {kz=0 plane) at the reduced Cg band hlling 
by 0.15 (hole doping): (a) La 2 Gu 04 (LDA), (b) La 2 Gu 04 (QSGW), (c) HgBa 2 Gu 04 (LDA), 
and (d) HgBa 2 Gu 04 (QSGW). The color represents the amount of da. 2 _y 2 and 
character. 


FIG. m (a) The three model parameters for the comparison to Ref. l32|. Note that AE, 


is denoted by ADs -|- 5Dt in Ref. l32|. (b) The difference of the calculated model parameters 
(Ecai) from the experiments (F^exp)- 


FIG. O The values of AEg^ (or ADs -|- 5Dt in the notation of Ref 32) estimated by LDA 


(blue squares), QSGW (red triangles), and RIXS data (green circles). 


FIG. ini The calculated AEe^ as a function of epitaxial strain. The ao=3.782 A is the 
experimental value for bulk La 2 Gu 04 . The compressive and tensile strain are simulated with 
a = 3.755 and 3.905 A considering the substrate of SrLaA 104 and SrTiOs, respectively j^ . 
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